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^ ' Abstract 

^ , We present general prescriptions for the asymptotic expansion of massive multi- 

Q \ loop Feynman integrals near threshold. As in the case of previously known pre- 

' scriptions for various limits of momenta and masses, the terms of the threshold 

expansion are associated with subgraphs of a given graph and are explicitly written 
^ \ through Taylor expansions of the corresponding integrands in certain sets of pa- 

\^ ' rameters. They are manifestly homogeneous in the threshold expansion parameter, 

so that the calculation of the given Feynman integral near the threshold reduces 
to the calculation of integrals of a much simpler type. The general method is 
illustrated by two-loop two-point and three-point diagrams. We discuss the use 
of the threshold expansion for problems of physical interest, such as the next-to- 
next-to-leading order heavy quark production cross sections close to threshold and 
matching calculations and power counting in non-relativistic effective theories. 



1. Introduction 



Many interesting processes in particle physics, in particular those in which heavy quarks 
participate, involve more than one mass scale. Such processes are notoriously difficult 
to calculate in perturbation theory beyond the one-loop level. To proceed one has to 
resort to approximations, either numerical or analytical. Among the latter, asymptotic 
expansions in certain ratios of mass scales appear most promising, because the analytic 
complexity is substituted by the algebraic complexity associated with obtaining a large 
number of terms in the expansion, which, however, can be delegated to a computer. 
To achieve this goal, the integrals that appear in the calculation of any given term in 
the expansion should (of course) be simpler than the original Feynman diagram. In 
particular, this means that the expansion should be manifestly homogeneous, that is, 
every integral that appears in the construction should contribute only to a single power 
in the expansion parameter. 

Let us consider, for the purpose of discussion, a quantity that depends on a single 
kinematic invariant and a particle mass m, such as the two-point functions of heavy 
quark currents, or the production cross section for a pair of heavy quarks. There exist 
general explicit prescriptions |jl|J^J^ (see for brief reviews) to obtain the asymptotic 
expansion of these (and other) quantities as — >■ oo or m — oo. In this paper we propose 
a prescription to obtain the asymptotic expansion close to threshold, that is, as — >■ 
Am?. This limit has not been explored systematically yet, although it is of considerable 
interest for a field theoretical description of non-relativistic systems. We illustrate our 
method through examples in Sects. 2 and 4. In these examples we consider only scalar 
propagators. The extension to fermion and gauge field propagators complicates the 
numerators of the integrals, but is straightforward methodically. The general structure 
of the expansion, which we formalize in Sect. 3, follows the same strategy used to obtain 
the the expansions in limits of large/small momenta and masses, although the threshold 
expansion is quite different, in particular because we have to deal with three different 
scales near threshold. The method obviously generalizes to the threshold production of 
two unequal-mass particles and in fact to any particle threshold, when some massive 
particles are slow, although we do not treat these cases in this paper. Physically, the 
mathematical problem of constructing an asymptotic expansion is closely related to the 
notion of effective field theories. For example, the homogeneity property mentioned above 
translates into the property of manifest power counting for the effective Lagrangian. 
The threshold expansion provides some insight into how to construct a non-relativistic 
effective field theory within dimensional regularization, an issue that has received some 
attention recently. We discuss this point, together with other conclusions, in Sect. 5. 

Recently Tkachov has also discussed the possibility of performing expansions near 
threshold 0. However, his prescription has been applied in only to the discontinu- 
ity of the 1-loop 2-point diagram with two masses. The non-trivial interplay of several 
small scales that characterize the threshold region is seen only in more complicated loop 
integrals. 
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Figure 1: One-loop vertex integral. Solid (wavy) lines denote massive (massless) propagators. 

2. Heuristic motivation and examples 

A typical explicit formula for the asymptotic expansion of a given Feynman integral Fr 
(corresponding to a graph F) in a given limit looks like 

Fr~5]A^^Fr, (1) 

7 

where the sum extends over a certain subset of subgraphs of the graph, and the operators 
Aij perform Taylor expansions in the variables that are small in 7. To arrive at this 
result, one can use the following heuristic procedure: 1. Determine the large and small 
scales in the problem. 2. Introduce factorization scales and divide the loop integration 
domain into regions in which each loop momentum is considered to be of the order of 
one of the scales in the problem. 3. Perform, in every given region, a Taylor expansion 
in the parameters, which are small in the given region. 4. After expansion, ignore all 
factorization scales and integrate over the entire loop integration domain in every region. 

One can easily reproduce the general formulae and combinatorical structure of the 
large-mass and large-momentum expansion in this way. The non-trivial point to justify 
is 4., which also guarantees the homogeneity of the expansion formula. In order for 4. 
to be valid it is essential to use dimensional (or analytic) regularization for the Feyn- 
man integral F^, even if F^ is finite in four dimensions. Loosely speaking, 4. follows 
in dimensional regularization from the property that all integrals without scale vanish. 
For off-shell limits of Feynman diagrams, the above procedure can indeed be justified 
in terms of the R* [^] and i?~^-operations. In this section we use these heuristic rules 
to treat two one-loop and one two-loop example and demonstrate that in each case the 
result agrees with the expansion of the exact result. 

Example 1. As our first example we consider a one- loop vertex integral with two equal 
and one zero mass, which appears for example in the form factor for 7* — > Q{pi)Q{p2) 
with quarks Q of mass m, see Fig. |l|. In the following we will also refer to massive 
lines as quarks and massless lines as gluons in general. We take Pi=P2 = and 
define q = Pi + P2 and the relative momentum p = {pi — P2)/2. The threshold region is 
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characterized by 



y = m - — = p <t:q . (2) 

We also define y = y/q^ as the dimensionless parameter of the threshold expansion. 
The parameter y is related to a more standard parameter /5 = (1 — Am? / q^Y^"^ through 
[3 = \J—Ay. Above threshold y < 0. It is useful to choose a frame in which q = (go, 0), 
Pi = {PojP), P2 = {po, —p) and p = (0,p). In this frame, the two massive particles move 
slowly and can be considered as non-relativistic. The following three scales are relevant 
to the threshold kinematics: q = ~ 2m, the centre-of-mass energy, ^—y = \p\, the 
relative momentum of the quarks and y/q ^ (Pi,2)o ~ the (non-relativistic) energy of 
the quarks. Accordingly, the loop momentum can be either 

hard (h): ko q, k q, 

or small (sm): ko<y/y, k<^. (3) 

Because there are two small scales, when the loop momentum is small, it can be either 

potential (p): kQr^y/q, k ^/y, 

soft (s): /i;o ~ v^, k^ (4) 

or ultrasoft (us): k^ ~ y/q, k ~ y/q. 

The terminology implies a 'canonical' routing of the large external momentum q through 
the massive lines of the graph only. For example, in Fig. ^ one should assign the mo- 
mentum fc±g/2 to the massive lines. Other routings are possible, but scaling arguments 
become less transparent then. The distinction between potential, soft and ultrasoft loop 
momentum is made only after, in the small momentum region, the integration over the 
zero-component of the loop momentum has been carried out by picking up the residues 
of the poles in /cq. The potential region is then associated with quark propagator poles 
and fco -C because the massive particles are non-relativistic. The soft and ultrasoft 
regions arise from the gluon propagator poles. We shall see that in general (but not in 
this example) the characteristic momentum of on-shell massless particles in the small 
momentum region can be either ^/y ~ mv or y/q ~ mf^, where v is the relative ve- 
locity of the two external quarks. (When the momentum of a massless line is in the 
potential region, the massless particle is off-shell by an amount of order y.) With these 
preliminaries the scalar integral of Fig. |l] is given by 

J _ / /.X 

^ J {k^ + q-k-y){k^-q-k-y){k-pY' ^' 

using the kinematic variables introduced above. The standard +iO-prescriptions are 
implicitly understood in the propagators and the integration measure is defined by [dk] = 
^e^E ^d^i (j^^d/2-\^ with d = A — 2e and 7^ = 0.577216 . . .. The renormalization scale of 
dimensional regularization is set to 1. 
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When the loop momentum k is hard, the integrand is expanded in y and p. (Because 
p ■ q = 0, terms odd in p vanish and the expansion produces only powers of y.) The 
leading term is 

= f m = (^TU) ^ (6) 

^ J k^{k^ + q-k){k^-q-k) \qy \2jl + 2e' 

Higher-order terms in y are calculated easily. The integral is evaluated most directly by 
reducing the number of propagators to two through partial fractioning. If one does not 
use partial fractions, the calculation of the integral with Feynman parameters leads to 
singularities inside the Feynman parameter integration domain when d = 4, which indi- 
cate the presence of the Coulomb singularity in /i as ?/ — > 0. This singularity is regulated 
by dimensional regularization as conventional ultraviolet and infrared singularities are, 
and /{* is well-defined. 

The contribution from the hard region corresponds to the 'naive' Taylor expansion of 
the integrand in y and p. However, the integral receives an important (in fact dominant) 
contribution from the small loop-momentum region, where (as we will verify shortly) the 
spatial loop momentum is of order ^/y. Inspecting the quark propagator, we see that it 
can now be expanded in k^: 

_J ^ ^ ^ I feo) (7) 

ki-k'^ + qoko - y „=o (-^^ + q^k^ - yY+^ 

The leading contribution to Ji from the small momentum region is then 

ir = I • (8) 

J (-F + go^o - y){.-k'^ - qoh - y){kl - {k-pf) 

We now perform the integration over /cq by closing the contour in the upper complex 
/co-plane. 

Let us consider first the contribution from the pole of the massive ('quark') propagator 
at fco = —{k"^ + y)/qQ + «0. Since \k\ ^ g by assumption, we have k^ -C \k\. As a 
consequence we might have expanded the massless propagator in (§) in k^ before picking 
up the residue from this pole. After fco-integration one obtains a two-point function with 
one mass y m d — 1 dimensions with the result 

lP = tR e^.. 1 = ...E JLL v^r(6+l/2) 

' q J T^'^'^-^ {P + y){k - pf V¥y 2e 

Because the integrand contains only one scale, the loop momentum is dominated by 
k ~ ^yy. This justifies the scaling given in (|^) for potential loop momenta. Note 
that higher-order terms in the expansion in k^ yield zero, because after integration over 
fco positive powers of ko result in massless tadpole integrals. Note also that we could 
have obtained (^ by expanding all terms that are small in the potential region before 
integration over ko- The contribution to the threshold expansion from the potential 
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region can then easily be estimated, because by construction all terms in the denominator 
have the same scaling in y. Taking into account that [dk] ~ y^^"^ /q for potential fc, we 
obtain If ~ \j \/q^y as expected for the Coulomb singularity and born out by (^). 

Now consider the contribution from the pole of the massless propagator in located 
at /co = —\k\ + iO. (We have shifted k ^ k + p. If one does not perform this shift, the 
combination k + p should be considered ultrasoft, which implies a cancellation between 
k and p, so that the scaling rules can not be applied to in a straightforward way.) The 
result is 

js/us ^ ^e-iE f ^'^ I QQ^ 

1 J ndn-^\k\[qQ\k\ + (k2^2p-k)][qo\k\-{k^ + 2p-k)]' 

Expanding in the small terms k"^ + 2p- k, the integral becomes a tadpole and is zero to 
all orders in the expansion. Hence, there is no contribution from the gluon pole. At first 
sight this looks incorrect, because the integral (0) is clearly non-zero. However, the non- 
zero contribution can come only from the region where k"^ + 2p- k is comparable to qo\k\, 
which requires k ^ q. This contribution is already included in the hard contribution 
above and (|TD|) must be expanded to avoid double-counting. Because the integral /^/""^ 
vanishes, we can not decide whether the gluon was soft or ultrasoft. 

There is a useful short-cut to arrive at this result. We can perform the approximations 
appropriate to the soft or ultrasoft region before integrating over ko to obtain, to leading 
order in this region, 

js/us ^ f [dk\ 

' J [qoko + tO][-qoko + tO]k^- ^ ^ 

This integral is ill-defined, because the poles at fco = pinch the integration contour. 
However, since the quark propagator poles have already been taken into account through 
the potential region, the previous integral should be understood as the contribution from 
the gluon pole only. Then one immediately arrives at a vanishing tadpole integral. Since 
by expanding all small quantities in the denominator, all remaining terms have the same 
scaling in y/y, we have homogeneity (manifest power counting) also for the soft and 
ultrasoft regions. 

One can now verify that the threshold expansion reproduces the expansion of the 
exact one- loop result: 

r(e) gr(l + e + n) 1 + 2e (-4y)"] 



2(l + 2e);^o r(l + e) 1 + 2e + 2n 



nl 



Note that in this example, contrary to those we discuss later, the Taylor expansions of 
the integrands do not generate additional poles in e. The poles in e in both the hard and 
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jh-p 




Figure 2: Diagrammatic representation of the threshold expansion for l2- Solid lines on the 
left-hand side denote massive propagators, wavy lines denote massless propagators. On the 
right-hand side a solid line with 'it' denotes (hard) 'on-shell' propagators of form l/{k'^ ± q.k) 
for line momentum kzizq/2, a wavy line stands for a hard massless line. The dashed line denotes 
a potential massive line and the dotted line a potential massless line. The zigzagged line is an 
ultrasoft massless line. 



the potential contribution come from the infrared pole which is present in the original 
integral Ji. 

Example 2. We now consider the two-loop 2-point master integral. This is the 
first case, relevant to the two-loop corrections to heavy-quark currents, where the hard, 
potential and ultrasoft regions contribute in a non-trivial way to the threshold expansion. 
The exact result for arbitrary y serves as a check of the expansion method. The diagram, 
shown in Fig. p|, is given by 



[dk][dl] 



{k^ + q-k- y){k^ -q-k- y){P + q-l- y){P -q-l- y){k - ly 



(13) 



and the lines in Fig. ^ are enumerated according to the order in which they appear in 
the denominator in (|T^). Expansions of this diagram in various limits of momenta and 
masses have been considered in [0. The 'canonical' routing of the external momentum 
for this diagram assigns ±g/2 to each massive line. With this assignment, the loop 
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momenta satisfy the scaling rules (0) and above. For the threshold expansion of any 
two-loop diagram, one considers the following nine loop momentum regions: h-h: both 
loop momenta are hard; h-p (s, us): one loop momentum hard, the other potential 
(soft, ultrasoft); p-p: both loop momenta potential; p-s (us): one loop momentum 
potential, the other soft (ultrasoft); s-s: both loop momenta soft; s-us: one loop momenta 
soft, the other ultrasoft; us-us: both loop momenta ultrasoft. In general, to obtain the 
contributions from all possible subgraphs it is necessary to consider different assignments 
of loop momenta than the one chosen in (0). As we discuss shortly, some of these regions 
do not contribute to the case at hand. The structure of the non-vanishing terms in the 
expansion is shown in Fig. 2. Note that I2 is finite in four dimensions. 

When both loop momenta are of order g, we can expand the integrand in y. After 
partial fractioning, the resulting integrals are of type 

J , ^ /■ [dk] [dl] 

J±^ai,...,a5j -y ^_^2Yi[-l2Y2[-[k-lYY^[-{k^ + q-k)Y^[-{P±q-l)Y^' ^ ' 

The integrals of type J+ can be reduced to Gamma-functions |^ through recurrence 
relations derived from integration by parts j^. The integrals of type J_ can be expressed 
in Gamma-functions and the integrals J_(0, 0, 03, 04, 05). These integrals can be reduced 



to J_(0, 0, 1, 1, 1), using a simplified version of the results of |jTO[- Then J_(0, 0, 1, 1, 1) 
(or a more convenient input integral) is calculated explicitly by Feynman parameters in 
an expansion in e. We obtain 



q-'lh-h = Q - 2 lng2 + 6 In 2j + 21C(3) - 4(8 + Svr'^) y + 0{f). (15) 

The h-h region is already the most difficult one. When one loop momentum is 
hard, and the other small, one of the three subgraphs 71 = {1, 2, 5}, 72 = {3, 4, 5} and 
73 = {1,2,3,4} can be hard. (The numbers refer to those in Fig. 2.) The lines in the 
hard subgraph 7^ are expanded in y and I (if k is hard) and the quark propagators not 
in 7j are expanded in the zero-components squared of momenta as in (^. The resulting 
contribution has the factorized form 

FTH^oTyF,^^ (16) 

i.e. the one-loop integral TyF^^, obtained by Taylor expansion of F^. in y and /, generates 
a local vertex, which is inserted into the reduced diagram r/7j obtained from F by 
shrinking all lines of 7^ to a point and expansion of the remaining massive propagators 
in the zero- components squared of the momenta. 

For our example, one obtains a vanishing integral for r/73 and equal contributions 
from 7i and 72. In the latter case, the non- vanishing contribution arises if the loop 
momentum of r/7i_2 is potential, but one obtains zero if it is soft or ultrasoft. The 
leading contribution from the h-p region is then (see Fig. 2) 

jh-p ^ 2 / ^3 / H (17) 
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where we expanded all terms that are small in the h-p region. It follows from power- 
counting that this integral contributes at order \fy. Keeping also the next term in the 
h-p region, we obtain 

g^j^^-p = Stt - lng2 _ Iny) 7^ - ^ fi - lng2 - Iny + f''' + O^f''). (18) 
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7i 



When / is soft or ultrasoft, the integrand is almost the same as for (|T^, except that / 
is expanded for ultrasoft / and P + y is expanded for soft I. In both cases, the resulting 
integrals vanish as already mentioned. 

When both loop momenta are small, we can use (|^) on all massive propagators 
and then take the poles in and /q- The resulting three-dimensional integrals are 
then divided into regions according to (|). As illustrated in example 1, one can also 
perform the appropriate approximations before integration over the zero components of 
loop momenta provided certain pinching quark poles in the soft or ultrasoft region are 
ignored as above. 

When both loop momenta are potential, all propagators are expanded in their zero- 
components squared. After picking up the residues from the remaining quark poles in 
/cq and Zo; the resulting d — 1-dimensional 2-loop integral is 

1 r d'^-^h d'^^^T 



^ 7rr(e + l/2)r(e-l/2) 
g2 2e 



,2 e''^^ ' ' ^ (19) 



as shown in Fig. 2. The next term and all higher order terms in the expansion in the 
p-p region vanish for the same reason as in example 1, hence 



q^lP-P = TT^ i^~-2 + A\n2 + 2\nyj . (20) 

We next consider the potential-ultrasoft region. Again we should consider all three 
cases, where one of the subgraphs 7j is potential and the remaining lines are ultrasoft. 
The only case that does not lead to vanishing scaleless integrals is potential 73, in which 
case the gluon line is ultrasoft. With k potential and / ultrasoft, the leading term is 



jp—us 

-'2 



[dk][dl] 

^-qo{ko-lo/2)~y) 

(21) 



(-P + q,{ko - lo/2) - y){-k^ - q,{k, - k/2) - y) 

1 



(-P + q,{k, + k/2) - y)i-k^ - qoiko + k/2) - y){ll - P) 

where the last (massless) propagator is now not expanded in Iq. (We have chosen a more 
symmetric loop momentum routing than in (|T3|).) Since [dk] ~ y^^"^ /q and [dl] ~ y^/q^-, 
power-counting tells us that this contribution scales as ^Jy. The integral is calculated by 
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examining the poles in ko and Iq and closing the integration contours in ko, /q such that 
the number of terms is minimized. One then arrives at 

^ ^ g2e7£; [ i ^22) 

2 ^2 J Txdl^-^ ^ J2 ^2{P + y) + qo\l\)' ^ ' 

Note that although \k\ ^ |/| in the region under consideration, the integral does not 
factorize in a manner comparable to ([T6|) , because of the presence of go ^ \k\ in the 
integrand. This can be phrased as the statement that integrating out the soft momenta 
does not yield a local interaction vertex with respect to the ultrasoft scale. This is in 
contrast to integrating out the hard momentum, which does result in a local interaction 
with respect to the small scales, see (]T6|), where TyF^- is polynomial in its (small) external 
momenta. Despite this fact, the above integral is effectively a succession of one- loop 
integrals. The result, keeping again the next term in the expansion in this region, reads 



,2; 



4^3/2^ 



g-/f-^^ = 87r(^---8 + 101n2-lng2 + 31nyj \^^y-Ji—)j ^ o{f'^). (23) 

A technical comment on the calculation of higher-order terms in y in the potential- 
ultrasoft region is in order. Expansion in the zero-components of the line momenta 
squared results in many terms, for each of which separately the integral over the semi- 
circle at infinity does not vanish upon closing the contours in the complex plane. As a 
consequence one has to take care of closing the contour in the same half plane for all 
terms, since the contribution from the semi-circle at infinity adds to zero only in the 
sum. 

Finally, when one momentum is potential and the other soft, we find scaleless in- 
tegrals. Likewise when both momenta are either soft or ultrasoft, one obtains only 
vanishing integrals. This concludes the list of all relevant momentum regions. 

In general, the h-h and p-p region contribute only to even powers of i/y and the 
h-p and p-us region contribute only to odd powers. The contribution from each region 
separately contains poles as e — >■ 0. The poles cancel between the h-h and p-p region and 
the h-p and p-us region, leaving logarithms of y. The terms computed above combine to 
the finite threshold expansion 

q^h = 27r2 ln(32y) + 21C(3) + 167r( ln(32y) - 4) - 4 (8 + Sn"^) y 

Wy)-f)f'^ + .... (24) 



The integral I2 is known exactly for any y. Using the result of |Tl|], we have 

q^h = F{l) + F{z^)-2F{z), (25) 

where 

F{z) = 6113(2) -41n;zLi2(-2) -\n^{z) \n{l-z), (26) 

1 - iy^ 



1 + 
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Pi 



P2 



P3 



k+pi-p3 



Pi 



Figure 3: The box graph. Sohd hnes massive propagators, wavy hnes massless. 

Taking care of correctly continuing the logarithms in to the second sheet when 

y < 1/4, the expansion of the exact result reproduces We have verified that this 

agreement persists to higher orders in y. 

Example 3. In the previous two examples, soft (on-shell) massless lines (gluons) did 
not contribute to the expansion. To obtain a non-vanishing result from this region, at 
least two gluons have to be exchanged between the massive lines (quarks). Consider 
a one-loop contribution to 2 ^ 2 scattering at small relative momentum. Fig. ^, and 
define pi/2 = q/2±p, P3/4 = q/2± p', where q + Pi + P2 = Ps + Pi, P = (pi - P2)/2 and 
p' = (p3 — p4)/2. All external lines are on-shell, pf = m^, and with y = — as 
before, we have = p'^ = y. We also introduce t = {p' — p)^ and choose a frame where 
q = (go, 0), p = (0,p) and p' = {0,p'). Then the threshold region is characterized by 

t^y<^q^. (28) 

The integral represented by Fig. || is 

(291 

{{k + pY + q-k-y){{k + pY -q-k-y)k'^{k + p-p'y' ^ ' 

The threshold expansion is obtained by letting the loop momentum be hard, potential, 
soft or ultrasoft and by expanding the integrand in all quantities that are small in the 
given region. In the following, we collect all terms up to order y^ jq^. (Recall that the 
tree diagram for 2 — > 2 scattering is 1/t ~ l/l/-) 

When the loop momentum is hard, the integrand is expanded in y, p and p' and one 
obtains 

When the loop momentum is potential, all four propagators are expanded in k\. This 
corresponds to picking up the quark propagator pole in the integration over /cq. Since 
the leading contribution scales as |/~^^^, subleading terms have to be included to achieve 
0{y^) accuracy in the threshold expansion. Defining t = t/q"^, the result is 

iq'r-I's = J^{\- H-i)) + Oif/M'^'). (31) 
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The subleading term actually vanishes in four dimensions. 

Let the loop momentum be ultrasoft. It is not possible to route the external momen- 
tum p or p' such that both gluon lines are ultrasoft simultaneously. Thus we get 

JUS _ 1 f [dk] 



'3 



1 / [^^1 = (32) 

t J [qoko + m-qoko + iOW ' ^ ^ 



where the pinch at ko — must be ignored, because it is taken into account by the 
potential region. 

However, there arises a non-zero contribution from the gluon poles, when the loop 
momentum is soft. In this case, to leading order in the soft approximation, the integral 
simplifies to 

r = i / [^^] (33) 

^ q^ J kik^{k + p-p'y ^ ^ 

Closing the /co-contour, we obtain 

r- ' [^""^ ' ( ' ' ] (34) 

3 q2 J p _ ^ ^ _ y^j^^s \%J^^_^>\Z) 

The singularity at P = {k -\- p — p'Y in the first factor is cancelled only in the sum 
of the two terms. Thus we cannot shift the momentum in the second term alone to 
make it equal to the first. We can calculate both terms separately, if we introduce a 
-|-iO-prescription in the first factor. After this we can shift the loop momentum and 
obtain 

^~ q^ J Ti'^l'^-^ (k'2f/2\-2k-{p-p')+t + iQ -2k-{p-p')+t-io) ' 

Now the integrals can be calculated by standard methods, taking care of the different iO- 
prescriptions in the two terms. Again keeping the first subleading term to the threshold 
expansion from this region, we find 

The contribution from the soft region is rather peculiar from the conceptual point of 
view. Both gluons in the box subgraph have energies of order ^/y, parametrically larger 
than the characteristic energy scale y/q for the non-relativistic quarks. Hence one can 
imagine the exchange of the two gluons as taking place on a much shorter time scale than 
the one that is characteristic for a non-relativistic system. The intermediate quark pair 
becomes off-shell by an amount \/q^ through the interaction with soft gluons, while the 
off-shellness of the quarks remains of order y in the potential and ultrasoft regions. The 
soft region thus can be considered intermediate between the hard region (off-shellness of 
order q'^) and the two other small momentum regions. However, while a hard subgraph 
can be interpreted as a local interaction vertex, the soft box subgraph, when it appears 
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as a subgraph in a 2-loop diagram gives rise to a vertex which is non-local in its external 
three-momenta but local in the zero-components of the external momenta as we discuss 
in more detail in Sects. 3 and 4. We also understand why the soft region did not appear 
in the previous two examples. Since an external quark pair or a quark pair that couples 
to the external current is off-shell by an amount of order y, after a quark interacts with 
a soft gluon which puts it off-shell by an amount of order y/q'^y, it has to interact with 
at least a second soft gluon to reduce its off-shellness to y again. 

The threshold expansion can be compared with the very simple exact (up to terms 
that vanish as e — > 0) result for the box diagram ||12|| : 



(9 



2\2+e : 



tVy 



TT 



— arctan J Ay 



1 



- In(-t) 



4 

i 



16y 
3t 



In(-t) 



(37) 



The expansion in y and i agrees with the sum of the terms computed above. Note 
the cancellation between the hard region and the constant from the soft region, which 
has to repeat itself in higher orders due to the factorized t-dependence of the exact result. 



3. General structure of the threshold expansion 

In the previous section we have been rather explicit and worked out the threshold 
expansion on simple examples. In this section we abstract from these examples general 
prescriptions for the threshold expansion {y = iv? — g^/4 and other small parameters 
such as i in example 3 tend to zero at fixed rn?) of an arbitrary diagram which contains 
two paths of lines with 'slow' particles with mass m ('quarks') and some massless par- 
ticles ('gluons'). (Internal loops of massive particles can be included straightforwardly: 
If all momenta that connect to the internal heavy quark loop are small, the integrand 
is expanded in external momenta over the mass and reduces to a series of local opera- 
tors analogous to the Euler-Heisenberg effective Lagrangian. If large momentum flows 
through part of the graph, the corresponding lines are not expanded. With this in mind, 
the following discussion is adapted to diagrams without internal heavy quark loops.) 
The two paths of massive lines are assigned q/2 and —q/2 of the large external momen- 
tum q. The two paths can be either disconnected as in the box graph example or they 
can be joined together by a hard external vertex that ejects large momentum q into the 
graph, such as in the two-point and three-point functions discussed in Sect. 2. We also 
assume a frame such that q = {qQ,0). The loss of explicit covariance is intrinsic to the 
threshold problem and, while an explicitly covariant formulation is possible, it would be 
exceedingly cumbersome. 

We arrive at our general prescription in three steps, which correspond to dealing 
with the hierarchy of scales sequentially. To arrive at the first form of the expansion, 
we consider any loop momentum to be either large or small. However, at this step, we 
do not identify small momenta as potential, soft or ultrasoft. In particular, when a loop 
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momentum that flows through a quark hue is small, the zero-component squared of the 
loop momentum in the quark propagator is small, see (|^. Thus, at the first step, we 
arrive at the following expansion 

^r = E%^r/,or,.,,F„ (38) 

7 

where the sum is in one-particle-irreducible (IPI) subgraphs of the given graph F. The 
operator T^i^y in the second factor performs the Taylor expansion of the Feynman integral 
in y and the loop momenta k'^ of F which are external with respect to 7. (These are 
the loop momenta of the reduced graph F/7, which is obtained by shrinking all lines 
of 7 to a point.) The operator 7^2^ in the first factor performs the Taylor expansion of 
the Feynman integral -Fr/7 associated with the reduced graph in the squares of the zero- 
components of the loop momenta that flow through the quark lines of the given paths. 
The above equation has the combinatorical structure of a factorized expression. It is 
easy to see that (|38[) is equivalent to constructing a non-relativistic effective Lagrangian 
TB| . The factors Tk-i^yF^ can be associated with the insertion of local operators into 



the diagram and they account for all hard loop momenta. The factors T^2^Fy/^ can be 
associated with diagrams in the effective theory and the hard scale q is not present as a 
scale for loop momenta in this part any more. Note that the asymptotic expansions for 



certain limits of on-shell integrals considered in can be recast into a form similar to 
(p8|), with expansion in k"^ rather than k^ in the small momentum part. 

However, the expansion (^8]) is not homogeneous in the expansion parameter. The 
non-homogeneity arises from the small momentum parts 7^2^ Fp/^, which are still non- 
trivial series in the threshold expansion parameter (s). To obtain the expansion in a 
manifestly homogeneous form we have to specify further the scale of the loop momenta 
of the reduced graphs F/7 in (|3^) . At this point we consider potential, soft and ultrasoft 
loop momenta with the scaling rules given in @j. Let F be a reduced graph that appears 
in (|38|). Since we have already disposed of all hard momentum regions, the 'next hardest' 
region is the region of soft loop momenta. Consider a soft subgraph 7^ of F, i.e. a 
IPI subgraph of F in which all loop momenta scale as /q" ~ ~ ^Jy and the 
external momenta k"'" of 7^ are either potential or ultrasoft. Then 7^ is expanded in all 
external momenta, which are ultrasoft, and reduces to a local interaction with respect 
to the ultrasoft scale. For a potential external momentum, we expand only in its zero 
component, because the spatial component of a potential momenta is of the same order 
as the spatial component of a soft momentum. As a consequence, the soft subgraph is a 
temporally local vertex with respect to the potential scale, but it is spatially non-local, 
because the dependence on the three-momenta of the potential external momenta is not 



polynomial. Therefore, defining Ff. = T^.2^Ff. as it appears in (p8D, we arrive at 

= E^f/,.°^e.,^-o^7., (39) 

7s 

where the sum is again in all IPI subgraphs. In addition to the Taylor expansions 
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indicated explicitly, the massive propagators in 7^ should be expanded according to 



and pinched poles are simply to be ignored, see the discussion in Sect. 2. Notice that 
these are static propagators rather than non-relativistic propagators. After inserting 
(^) in (p^) we have achieved factorization of both hard and soft momenta. The factors 
"^kZl ,k''%-^'ys ^^'^ be interpreted as an instantaneous interaction vertex, and thus as con- 
tribution to the heavy quark potential. On the other hand, since Tj^is^^j^is^F^^ is spatially 

non-local, the graph r/7s is interpreted as deleting all lines of 7^ in F and replacing them 
by 7^7s ;i,7s F^^. Let us give an example. The soft part of the box graph of Sect. 2 has 
four potential external quark lines and we have 

c 1 

r^x^F,^. - ^ F(^^7?F ^^^^ 

with a numerical constant c. In coordinate space, this gives rise to a non-local four- 
fermion operator. Due to the structure of (|39|) , the factorization of soft regions can 
also be implemented as an effective Lagrangian, although it is non-local in space. In the 
planar two-loop three-point integral which we treat in Sect. 4, we will see how such a non- 
local operator inserted into a reduced diagram contributes to the threshold expansion 
of two-loop diagrams. In Sect. 4 we will also meet another non-local vertex for the 
interaction of two potential quark lines with one potential gluon line. One can use power 
counting for the soft region to determine all non-local operators which contribute to a 
given order in the threshold expansion. 



The expansion obtained from combining (pSD and (|39D is still not homogeneous in the 



threshold expansion parameter. The inhomogeneity arises from F(r/^y^^, where hard and 
soft lines are deleted, but potential and ultrasoft loop momenta are still not separated. 
First, we note that quark lines can never be ultrasoft. Technically, this is so, because, 
since only ko can combine with a large qo, spatial tadpoles always result, if one tries 
to make quark lines ultrasoft. As a consequence, one should consider all contributions 
in which all quark lines are potential and a (possibly empty) collection of (possibly 
disconnected) gluon lines is ultrasoft. When an ultrasoft gluon line with momentum / 
connects to a quark line with loop momentum k — 1/2 for the incoming and k + 1/2 for 
the outgoing quark line, we can expand the quark-gluon vertex and quark propagator 
in //A; ~ ^/y. This corresponds to a multipole expansion, as the wavelength of ultrasoft 
gluons is large as compared to separation of the quark-antiquark system in a coordinate 
space picture. After multipole expansion, all scales are separated and we arrive at a fully 
homogeneous threshold expansion. 

We have arrived at this prescription in a heuristic and descriptive way. We are 
presently not able to give a mathematical proof of our prescriptions similar to what 
is available for off-shell limits and leave a mathematical treatment to future pub- 
lications. To complete such a proof for the threshold expansion one needs to use the 
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a-representation and a subsequent decomposition of the integration domain into appro- 
priate sectors which should be introduced, in the language of the a-parameters, in a 
way analogous to our decomposition into regions of large, potential, soft and ultrasoft 
momenta. In fact, a proof of the asymptotic expansions for off-shell limits reduces 
to the analysis of the asymptotic behaviour of the remainder, which has the structure 
of the /^-operation, i.e. renormalization at the diagrammatic level. Therefore this proof 
has much in common with standard proofs of the finiteness of the renormalized Feynman 
integrals. 

According to the definition of the asymptotic expansion, its remainder possesses a 
given asymptotic behaviour (a sufficiently fast decrease) once one keeps a sufficiently 
large number of the terms of the expansion. As in the case of the limits for which 
explicit general formulae were known up to now [p],p|,p|,p!4| , it is possible to characterize 
the remainder of the threshold expansion in terms of the /^-operation. Since the product 
of any two (or more) subtraction operators that are present in the prescriptions for the 
asymptotic expansions is zero, we can rewrite the sum that in the right-hand side of the 
expansion (|ip as 

Y^M^ = l-n (42) 

7 

where the operation TZ has the structure of the usual (ultraviolet) /2-operation and is 
given by the forest formula 

^ = E n (--^7) • (43) 

For off-shell limits, the operator Ai^ performs Taylor expansion of the Feynman integral 
in its small masses and external momenta and inserts the resulting polynomial in the 
reduced diagram F^/^. In our case, the Taylor expansions appropriate to a given region 
should be implied. In ( ^31 ) the sum is in forests (i.e. subsets of non-overlapping sub- 
graphs) composed of subgraphs that are involved in the prescriptions for the asymptotic 
expansion in the given limit. Therefore, we can represent an initial Feynman integral as 

Fr = J2 -^7^r + . (44) 

7 

The remainder is constructed in such a way that it is finite (if the original diagram 
is finite) and possesses a desired estimate in the given limit. The order of the expansion 
is implicitly hidden in degrees of the subtractions operators involved. From the mathe- 
matical point of view, the subtractions performed by the operators corresponding to 
subgraphs 7 7^ F remove divergences generated by the naive expansion (i.e. for 7 = F) 
of the integrand in the expansion parameter. This naive expansion in turn removes 
divergences that could appear due to the subtractions in the subgraphs. 

Let us illustrate the above points through our examples. In example 1, we have the 
following expression for the remainder: 

. jm (1 - T/) (1 - T|) ^ ^ . , . ^)(,. , . , _ _ ■ (45) 
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Here the superscript denotes the order of the Taylor expansion. It turns out that 
the remainder possesses a necessary asymptotic behaviour, provided the degree of the 
subtraction operators is chosen sufficiently large, and does not have divergences that 
were not present from the very beginning. Consider first the region of small k. Then 
the operator Ty^ is dangerous, in the sense that it generates an infrared divergence that 
appears when y 0. However, the second factor (1 — 7^^) kills this divergence because 
it either provides additional factors of A;, or explicit powers of y (which effectively reduce 
the order of expansion in y in the first factor). Similarly, in the region of large k, the 
operator T/^ is dangerous because it generates positive powers of k. This time, the first 

factor (1 — 7^^) improves the convergence of the whole integral. 
Consider example 2. As a remainder we now have TZF^ where 

7^ = (1 - Mr)il - Mh-p - Mp-h)il - Mp-p - Mp-us) ■ (46) 

The operators M.h-p, M.p-h, -Mp^p and M.p_us have to remove infrared divergences 
generated by the action of the operator A^r = -Mh-h = %■ Consider, for example, the 
region where k is small while I is non-zero. Then one can pick up the bracket (l — Aih-p) 
which makes the term with A^r-^r convergent. In the region where both k and I are 
small, this is the operator Mp-p that provides convergence that could be spoiled by the 
operator A^r, and the operator A4p_us that kills the divergence at /c, Z — > caused by 
Mh-p and Mp^h- 

As in the case of known limits of momenta and masses, we have, in the right hand- 
side of (|45|) and (^6]), and ( pH]) and (|39|), an interplay between ultraviolet and infrared 
divergences which manifest themselves as poles in e. However, the naive expansion (in y) 
does not always generate additional poles — see, e.g., example 1. Still the naive Taylor 
expansion part taken alone does not provide a correct answer to the problem and should 
be accompanied by extra contributions. 



4. Two-loop vertex integrals 



In this section we discuss and present results on the threshold expansion of two-loop 
3-point functions, which, to our knowledge can not be obtained from expansion of an 
exact result. The three-point functions are also the non-trivial building blocks for many 
physical applications. Here we consider the abelian planar and non-abelian graphs in 
Fig. ^ and b. The non-planar graphs (Fig. ^) will be discussed together with first 
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The abelian planar graph. Define the one-loop subgraphs 71 = {1,3,5}, 72 = 
{2,4,5,6} and 73 = {1,2,3,4,6}. (The numbering of lines refers to Fig. ^.) The 
unexpanded integral is given by 



'PL 



[dk][dl] 



{k^ + q-k- y){k^ -q-k- y){P + q-l- y){P -q-l- y){k - Ifil - pf 

(47) 
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(a) (b) (c) 

Figure 4: Examples of two-loop three-point graphs. Solid (wavy) lines denote massive (mass- 
less) lines. 

where the kinematic variables are defined as in example 1 above. Inspecting all regions 
that contribute to the threshold expansion, we find only six of them non-vanishing. 
Together with their leading contribution (up to logarithms of y) to g^/pL, they are: 

h-h: 1 
h-p, with 7i hard: 

h-p, with 72 hard: (48) 
p-p: 1/y 

p-s, with 72 soft: l/\fil 
p-us, with 73 potential: l/yy. 



Let us emphasize that the scalings given here follow from power counting alone, because 
in each region the integrand is homogeneous in the expansion parameter. 

The calculation of the various contributions is straightforward in the sense that all 
non-trivial aspects of the calculation have been met in the examples discussed in detail 
in Sect. 2. For instance, after partial fractioning, the integrals needed to evaluate the h-h 
region reduce to J± of (|T^ (with additional scalar products in the numerator) and the 
complexity of integrals is essentially the same as for the 2-point function of example 2. 
The leading contribution close to threshold comes from the p-p region and corresponds 
to the double Coulomb exchange. Let us consider in more detail the fifth contribution, 
when the box subgraph is soft and loop momentum k that flows through lines 1, 3 and 5 
is potential. First, when we take the integral over the other loop momentum / first, the 
result is identical to the soft contribution to the box subgraph in example 3 of Sect. 2, 
except that the difference of external relative momenta p—p' is replaced by the difference 
of the external momenta of the box subgraph p + k—p = k. The momentum dependence 
of the box subgraph is 1/(A;^)^^^, which can be interpreted as a non-local vertex. Second, 
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using p ■ g = 0, it is easy to see that in every region the expansion parameter is y, except 
when one momentum is soft, in which case the expansion parameter is \/y. However, for 
the box subgraph the symmetry between the upper and lower quark hne ehminates all 
odd terms in -y/y and subleading contributions are suppressed by integer powers of y. 

The integral Jpl is infrared divergent. Collecting all terms up to order y° close to 
threshold, we obtain 
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(49) 



The calculation of higher-order terms can be automated on a computer. As a check of 
this result, we have calculated the exact double pole part of the integral. 
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and the first three terms in the expansion agree with the double pole part of 



(50) 



The non-ahelian graph. For the diagram of Fig. we define the one-loop subgraphs 
7i = {2, 4, 5}, 72 = {1, 3, 4, 6} and 73 = {1, 2, 3, 5, 6}. The unexpanded integral is given 
by 



'NA 



[dk][dl] 



{k^ + q-k- y){k^ -q-k- y){P + q-l- y){k - If^k-pYil - 



(51) 



Inspecting all regions that contribute to the threshold expansion, we find only four of 
them non- vanishing. Together with their leading contribution (up to logarithms of y) to 
g^/NA, they are: 

h-h: 1 

h-p, with 7i hard: l/-y/^ (52) 
p-s, with 7i soft: 1/y 

p-us, with 73 potential: l/\/^. 

The p-p region that scales as 1/^ does not contribute here, because the poles in Iq lie on 
one side of the real axis only. Hence, the /o-integral vanishes. The potential-soft region 
is algebraically most complicated. Since the soft vertex subgraph 71 is less symmetric 
than the box subgraph of the abelian planar diagram, the expansion in this region runs 
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in \fy and it is necessary to compute the first two subleading terms in the expansion 
in this region to obtain an accuracy 0[ip^ at threshold. The integral Jna is infrared 
divergent. Collecting all terms up to order \p close to threshold, we obtain 
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As a check of this result, we have calculated the exact double pole part of the integral, 

1 1 /TT 




TT / \ 

— — arctan y 4y j , (54) 



IGe^Vy 1 + 4y V2 

and the first two terms in the expansion agree with the double pole part of (l53|) . 

We close this section with some observations on the non-planar diagram. Fig. |c. 



postponing a detailed discussion to [15]. The non- vanishing regions are the same as for 
the non-abelian diagram, with the difference that the hard-potential contribution scales 
as a/^ in leading order and the potential-soft region as 1/ \fy as for the planar abelian 
diagram. (The p-p region vanishes again, because all poles in one of the zero- component 
integrations lie on one side of the real axis.) Any integral that contributes to the h-h 
region can be reduced to either J± or 



L±(ai, ... ,05) = 

\dk\\dl\ 



(55) 



With integration by parts we can reduce these integrals first to L+(0, 0, 03, 04, as) = 
L„(0, 0, 03, 04, as) and then to a single integral such as -L+(0, 0, 2, 2, 1). (The algorithm 
for L+ has been used before in P,p!6|, although no details have been given there.) Fur- 
thermore, the hard parts of all possible 2-loop 3-point graphs relevant to QCD can be 
reduced to J± or L± and, consequently, can be solved in the threshold expansion. 

5. Discussion 

The threshold expansion, which we proposed and illustrated in this paper, could be 
applied to two sets of problems. The first set consists of problems where it is desirable to 
have as many terms in the expansion as can possibly be calculated. The second consists 
of problems where the nature of the problem specifies that only a limited number of 
terms are needed. 

The cross section for the production of heavy quarks in e^e~ annihilation o"e+e- [(^ ■, ttl) 
at intermediate energies \/q^ is an important example of the first kind. Presently avail- 



able techniques do not permit an analytic calculation of o"e+e- {q'^, fn) at order a^. In [17 
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the large-momentum and large-mass expansion together with convergence accelerating 
methods have been used to approximate cXe+e- (o'^, ^) to high accuracy from expansions 
around = and = oo. It is exactly the threshold point = Am? that is most crit- 
ical in this approach. The threshold expansion, when computed to high order, provides 
us with an intermediate point, the expansion around which can probably be smoothly 
joined with the large-momentum and large-mass expansion in order to obtain an accurate 
result for (Je+e-{q^i^) for any q^ except very close to threshold, where a resummation 
of all corrections of form Q;°'^'^(a<i/v^)" is necessary. (This resummation belongs to the 
second set of problems.) The analysis of the 2- loop 2-point function in Sect. 2 suggests 
that the threshold expansion for a^+^-{q'^,m) converges when G [0,8m^]. Therefore 
we expect that a reasonably accurate approximation to ae+e-iq'^yfn) could be obtained 
for < 6m^ starting from the threshold expansion. 

In this paper we have treated only the threshold expansion of loop integrals, but not 
of phase space integrals for real radiation. However, as in it seems advantageous to 



consider the three-loop vacuum polarisation of heavy quarks in the threshold expansion 
in order to avoid the separate calculation of virtual corrections and real radiation. Since 
(jg+e- (g^, "O^) requires only the imaginary part of the two-point function, one can take 
advantage of this fact from the beginning. For example, the region where all three 
loop momenta are hard, which would be very difficult to calculate, is analytic in y and 
therefore does not contribute to the imaginary part. Other contributions reduce to two- 
loop diagrams (such as those of Sect. 4) or simple three- loop diagrams. 

The second set of problems involves 'matching calculations'. In this case the quark 
mass m is a large scale and the observable in question should be factorized into a short- 
distance contribution from the scale m and a long-distance contribution which can be 
either perturbative or non-perturbative. Examples of this kind include bound state cal- 
culations in QED, quarkonium systems in QCD, but also the above-mentioned resum- 
mation of 'Coulomb-enhanced' corrections in heavy quark production close to threshold. 
In order to obtain the matching coefficient for a given operator, a suitable (typically 
on-shell) Green function has to be expanded close to threshold to a certain order, which 
is determined by the structure and dimension of the operator. As already mentioned in 
Sect. 3, the combinatorical structure of the separation of large (hard) and small momenta 
in the threshold expansion, which is manifest in ( pS]) for a given diagram F, translates 
into the statement of factorization 

G{m,pi) = ^Ci{m)Go,{Pi/m) (56) 

i 

for the Green function G to which F contributes. Here Oi is a local operator, inserted into 
G, and Cj is its short-distance coefficient. This construction is completely equivalent to 
the construction of a non-relativistic effective theory |]T3| in dimensional regularization. 



In particular ( pHD leads to the factorization prescription that for any diagram F its 
contribution to the short- distance coefficient(s) ('hard part') is given by selecting only 
7 = F in the sum over all subgraphs: 

^rihard = Vr- (57) 
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That is, the hard part is given by the 'naive' Taylor expansion of T. Note that this 
is not identical to dropping all singular terms at threshold. For example, the abelian 
planar 3-point integral in Sect. 4 contains terms proportional to both from the h-h 
and the p-p region. The contribution from the second region should not be considered 
as part of the short- distance coefficient. The non-relativistic approximation to the quark 
propagator in does not coincide exactly with the standard convention. Let the 
gluon line in the triangle graph of Fig. |l| have momentum k and the upper quark line 
have momentum k + pi. In the standard convention one would introduce an external 
momentum p = {po,p) with the non-relativistic on-shell condition po = p"^ /{2m). The 
non-relativistic propagator is then 

^ '-^ = I . (58) 

*"po + *o-^J^ -{k + p)-'-+q«h-y + l> 

where the difference to the approximation used in this paper is 5 = (go/2 — m) /cq = 
(\/m^ — y — m) /cq) which is smaller than all other terms in the denominator and can be 
expanded. Thus our small-momentum approximation is equivalent to a reparametriza- 
tion of NRQED/NRQCD, where the large part of the massive line momentum is taken 
to be q/2 rather than m. 

Once the hard momenta are factorized ('integrated out'), an interesting aspect of 
non-relativistic effective theory is that it does not lead to manifest power counting. This 
means that the matrix elements of operators are series in the small parameters, rather 
than being homogeneous, and their scaling can not be determined by counting a suitable 
dimensionful parameter. In the threshold expansion this is reflected by the fact that we 
had to divide the small momentum region further into potential, soft and ultrasoft (see 
(D) in order to arrive at a homogeneous expansion with manifest power counting. 

Recently the possibility to factorize the small-momentum scales mv and mv"^ (y^y 



and y/q in our notation) has been investigated in a series of papers [pr8| , p!9| , pO| , pl| 



Different rescalings were proposed in lTyj2^ that succeed in making power counting 



manifest in one region, but fail in others. Thus while the rescaling of ||i9[ is adequate 
for potential gluons (/cq ~ mv^, k ~ mv), the rescaling of |^ is tuned for ultrasoft 
gluons (fco ~ /c ~ mv'^) and the multipole expansion. Both rescalings are unified in the 
treatment of where different fields for (non-propagating, off-shell) potential gluons 
and (propagating) ultrasoft gluons are introduced. At least for diagrams, in which only 
one gluon (massless) line is present, this approach appears to be equivalent to ours, as 
can be seen by comparing the treatment of the small momentum contributions to the 
triangle graph in Sect. 2 with the discussion in I^Tj. Furthermore, in the third line of 
Fig. 1^ the first graph corresponds to the exchange of a potential gluon and the second 
to the exchange of a propagating gluon in the language of . On the other hand, the 



existence of soft gluons, i.e. propagating gluons with energy and momentum of order mv, 
has not been considered in [^,^,^ and does not seem to fit into the rescaling schemes 
proposed there. Yet another approach has been taken in WM, where a non-local effective 



Lagrangian is proposed that follows after integrating out soft gluons. (In practice, this 
has been done only for single gluon exchange, which gives rise to the standard Coulomb 
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potential with relativistic corrections. The discussion of the soft box subgraph in Sects. 2 
and 4 suggests that the non-local Lagrangian can be extended to include the box graph 
contribution.) However, it has not yet been shown how to treat potential gluons in the 



approach of [22 



Labelle has first discussed the separation of soft / potential and ultrasoft photons and 
the multipole expansion in the context of NRQED It seems to us that Labelle's 

treatment of NRQED graphs is most closely related to our treatment of the small mo- 
mentum regions in the threshold expansion. Technical differences arise in as much as he 
implies a cut-off regularization, while the simplicity of the threshold expansion discussed 
here is specific to dimensional regularization. In addition, |jl8[ starts with time-ordered 



perturbation theory and Coulomb gauge, while we emphasize that the power counting 
and approximations in a given region can be carried out before performing the inte- 
gration over the zero- components of the loop momenta and the method works in any 
gauge. (Since we have considered scalar propagators only, the question of gauge does not 
arise. In a gauge theory Coulomb gauge is very convenient for the threshold expansion, 
because propagating gluons which cause most complications couple to quarks through 
vertices that are suppressed by at least a factor of ^y.) In this paper, as in and 



in contrast to |fr9ip0tpi] , P2| , no attempt is made to write down an effective Lagrangian 
after integrating out soft or soft together with potential gluons. This may be considered 
as a drawback. On the other hand, as this Lagrangian would inevitably be non-local, its 
utility is not obvious to us, once manifest power counting can be achieved otherwise. 

After this conceptual postludium we should emphasize that the main focus of this 
paper is on the ability to perform calculations. The threshold expansion method de- 
scribed here allows us to treat 2-loop 3-point functions and opens the door to threshold 
problems at next-to-next-to-leading order. A first realistic application will be treated in 
a companion paper ||15 



Acknowledgements. We thank A. Czarnecki, A. Davydychev and K. Melnikov for 
helpful discussions and providing checks on some of the integrals required for the hard- 
hard region. We are also grateful to A.V. Kotikov and I.Z. Rothstein for useful discussions 
and to W. Beenakker and B. Tausk for correspondence regarding the exact box graph. 
V.S. thanks the CERN Theory Group and the Universities of Leiden and Mainz for 
hospitality during the course of this work. Both, V.S. and M.B. thank the Theory 
Group at the University of Karlsruhe for hospitality while this work was completed. 
Thanks to G. Buchalla and K. Chetyrkin for reading the manuscript. The work of V.S. 
has been supported by the Russian Foundation for Basic Research, project 96-01-00726, 
and by INTAS, grant 93-0744. 



22 



References 



[1] 

[2] 
[3] 
[4] 

[5] 
[6] 
[7 



[9 
[10 

[11 
[12 

[13: 

[14 

[is: 

[16 

[17: 



S.G. Gorishny, Nucl. Phys. B319 (1989) 633. 

K.G. Chetyrkin, Teor. Mat. Fiz. 75 (1988) 26; ibid. 76 (1988) 207; Max-Planck- 
Institute preprint MPI-PAE/PTh 13/91 (Munich, 1991). 

V.A. Smirnov, Commun. Math. Phys. 134 (1990) 109; V.A. Smirnov, Renormaliza- 
tion and asymptotic expansions (Birkhauser, Basel, 1991). 

V.A. Smirnov, Mod. Phys. Lett. AlO (1995) 1485; Talk given at 5th International 
Conference on Physics Beyond the Standard Model, Balholm, Norway, 29 Apr - 4 
May 1997 ||hep-ph/970842B| . 



F.V. Tkachov, INR-F5T-97-01 [|hep-ph/9703424 



K.G. Chetyrkin and V.A. Smirnov, Phys. Lett. 144B (1984) 419. 

D.J. Broadhurst, J. Fleischer and O.V. Tarasov, Z. Phys. C60 (1993) 287; A.L Davy- 
dychev, V.A. Smirnov and J.B. Tausk, Nucl. Phys. B410 (1993) 325; F.A. Berends, 
A.I. Davydychev, V.A. Smirnov and J.B. Tausk, Nucl. Phys. B439 (1995) 536; 

F. A. Berends, A.I. Davydychev and V.A. Smirnov, Nucl. Phys. B478 (1996) 59. 

N. Gray, D.J. Broadhurst, W. Grafe and K. Schilcher, Z. Phys. C48 (1990) 673. 
K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159. 
O.V. Tarasov, DESY-96-240 ||hep-ph/97033Tl . 

A. Djouadi, Nuovo Cim. lOOA (1988) 357; D.J. Broadhurst, Z. Phys. C47 (1990) 
115. 

W. Beenakker and A. Denner, Nucl. Phys. B338 (1990) 349. 

W.E. Caswell and G.P. Lepage, Phys. Lett. B167 (1986) 437; B.A. Thacker and 

G. P. Lepage, Phys. Rev. D43 (1991) 196; G.P. Lepage, L. Magnea, U. Magnea, 
C. Nakleh and K. Hornbostel, Phys. Rev. D46 (1992) 4052. 

V.A. Smirnov, Phys. Lett. B394 (1997) 205; A. Czarnecki and V.A. Smirnov, ibid. 
B394 (1997) 211; V.A. Smirnov, Phys. Lett. B404 (1997) 101. 

M. Beneke, A. Signer and V.A. Smirnov, in preparation. 

J. Fleischer and O.V. Tarasov, Comp. Phys. Commun. 71 (1992) 193. 

K.G. Chetyrkin, J.H. Kiihn and M. Steinhauser, MPI-PHT-97-029 



ph/9705254|| . 



[18] P. Labelle, MCGILL-96-33 ||hep-ph/960849T 



23 



[19] M. Luke and A.V. Manohar, Phys. Rev. D55 (1997) 4129. 

B. Grinstein and I.Z. Rothstein, UCSD-97-06 ||hep-ph/ 9703291 . 



[21] M. Luke and M.J. Savage, UTPT-97-12 ||hep-ph/97073T^] , to be published in Phys. 
Rev. D56. 

[22] A. Pineda and J. Soto, UB-ECM-PF-97-17 ||hep-ph/9707481^ , to be pubhshed in 
the proceedings of High-Energy Physics International Euroconference on Quantum 
Chromodynamics: QCD 97, Montpellier, France, 3-9 Jul 1997. 



24 



